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Abstract 



We report on a search for a supersymmetric B meson with mass between 
3.5 and 4.5 GeV/c 2 using 4.52 fb _1 of integrated luminosity produced at 
y/s = 10.52 GeV, just below the e + e~ — ► BB threshold, and collected with 
the CLEO detector. We find no evidence for a light scalar bottom quark. 
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There has recently been renewed interest in the possibility of a light scalar bottom 
quark [l], which ALEPH has searched for ||. It has been noted that such a squark can 
exist in certain regions of parameter space [[J. In addition, the production cross section 
of scalar quarks in e + e~ annihilations, well above threshold, is \ that of spin-| quarks of 

the same charge; thus bb production would contribute ^ unit to R, the ratio of hadronic 
cross section to fi + fi~ cross section, and such an increase cannot be ruled out by existing 
measurements ||. Indirect evidence from the measured B semileptonic branching fraction 
disfavors the existence of a light b ||. 

If b is the lightest supersymmetric particle (LSP), and if .R-parity is conserved, then b 
would be stable. If, instead, a scalar neutrino v is the LSP, then b will decay b — > cl~v 
and/or b — > ui~v. If i?-parity is violated, b will decay b — ► ci~ and/or b — > ul~ . We have 
searched for a light b that decays b — ► c£~z> and/or 6 — ► c^~. Such a particle would dress itself 
as a supersymmetric S meson. The dressed decays would be B — > DX£~i> and B — > DX£~ , 
where X represents possible additional conventional hadrons. 

The decays we search for, characterized by leptons and charmed mesons, have much in 
common with conventional B decays. We perform a direct search, avoiding the B background 
by using a data sample collected below the BB threshold, at y/s = 10.52 GeV. Our search 
covers the B mass range 3.5-4.5 GeV/c 2 . Because our search is near BB threshold, the 
production cross section cannot be predicted with great precision. We include the /3 3 thresh- 
old factor and ignore strong interaction effects to obtain a(e + e~ — > bb) = N c na 2 Q'ff3 3 /3s 

(derived from Eq. 35.12 of Ref. |J), where Qf, the charge of the b, is — |. 

Below BB threshold, the major source of leptons and charmed mesons is e + e~ — > cc, 
which predominantly produces two charged leptons when both charmed mesons decay 
semileptonically. To suppress this cc background, we search for events with two oppositely 
charged leptons as well as a fully reconstructed hadronic D or D* meson decay, where D^*> 
denotes either a or £)(*)+ meson. Other sources of leptons include kaons that decay 

in flight, photon conversions, and 7r° Dalitz decays. A search for wrong-sign DW-lepton 
combinations was also conducted, but the resultant upper limits were significantly weaker 
than for the Z}(*)-dilepton signature, and they are not discussed further in this paper. 

The data sample used in this analysis was produced in symmetric e + e~ collisions at 
the Cornell Electron Storage Ring (CESR) and collected with the CLEO detector in two 
configurations, known as CLEO II J7[ and CLEO II. V ||. It comprises 4.52 fb _1 of integrated 
luminosity produced at an e + e~ center-of-mass energy of y/s = 10.52 GeV, 35 MeV below 
the BB threshold. Analysis of an additional 2.24 fb _1 collected on the T(AS) resonance 
provides a sample of B mesons, which can mimic the behavior of B mesons, thus allowing us 
to verify our experimental technique. The response of the experimental apparatus is studied 
with a GEANT-based || simulation of the CLEO detector, where the simulated events are 
processed in a fashion similar to data. 

In CLEO II, the momenta of charged particles are measured with a tracking system 
consisting of a six-layer straw tube chamber, a ten-layer precision drift chamber, and a 51- 
layer main drift chamber, all operating inside a 1.5 T superconducting solenoid. The main 
drift chamber also provides a measurement of specific ionization energy loss, which is used 
for particle identification. For CLEO II. V, the six-layer straw tube chamber was replaced by 
a three-layer double-sided silicon vertex detector, and the gas in the main drift chamber was 
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changed from an argon-ethane to a helium-propane mixture. Photons are detected with a 
7800-crystal Csl electromagnetic calorimeter, which is also inside the solenoid. Proportional 
chambers placed at various depths within the steel return yoke of the magnet identify muons. 

Charged tracks are required to be well-measured and to satisfy track quality require- 
ments based on the average hit residual and the impact parameters in both the r-0 and r-z 
planes. Muon candidates must penetrate the steel absorber to a depth of at least five nuclear 
interaction lengths, which effectively places a lower bound on the muon momentum of 1.2 
GeV/c. Electrons are identified by a likelihood that includes the fraction of the particle's 
energy deposited in the calorimeter and the spatial distribution of the deposited energy. To 
reduce contamination from low-momentum hadrons, the electron momentum is required to 
be greater than 1.0 GeV/c. For electrons that are combined with another lepton satisfying 
the above criteria, the momentum requirement is lowered to 600 MeV/c. Electron pairs from 
photon conversions are rejected by requiring the dielectron invariant mass to be greater than 
200 MeV/c 2 . We rely on the detector simulation to estimate the remaining contribution 
of hadrons to the lepton sample. 7r° candidates are formed from pairs of photons with an 
invariant mass within 2.5 standard deviations of the known 7r° mass that are kinematically 
fitted with the mass constrained to the known 7r° mass. 

D mesons are reconstructed in the modes D° —>■ K~tt + and D + — > K~7t + tt + , where a sum 
over charge conjugate modes is implied. The daughters of the D candidates must undergo 
ionization energy loss consistent with the particle hypothesis at the level of three standard 
deviations. To maximize detection efficiency, no requirements are placed on the momentum 
of the D candidate. We reconstruct D* mesons in the modes D*° — > D°7T , D* + — > D tt + , 
and D* + — > D + tt°. When used in the D* modes, D candidates must have invariant masses 
less than 15 MeV/c 2 (about 2.5 standard deviations) from the known mass. 

Figure [l] shows the K~n + and K~tt + tt + invariant mass distributions of D° and D + 
candidates, respectively, associated with two oppositely charged leptons. Also shown are 
the normalizing distributions for D candidates paired with a single lepton of either charge. 
The analogous distributions for the D* modes, showing the D*-D mass difference, are given 
in Figure [2|. The D^^-dilepton distributions reveal a striking absence of signal. The 
yields are extracted from a fit of each histogram to a Gaussian distribution over a linear 
background for the D modes and a quadratic background for the D* modes. The means 
and widths of the Gaussians for the Z)(*)-dilepton fit are fixed to values determined from 
the normalizing distributions, where the Gaussian widths are typically 6.5-7.0 MeV/c 2 for 
M(D) and 0.5-1.0 MeV/c 2 for M(D*) - M(D). Table | lists the yields observed in data 
as well as the expected yields determined from a simulation of e + e~ — > qq events, where 
q G {u, d, s, c}, hereafter referred to as "continuum" events. No significant excess is observed 
in any mode. 

In the normalizing distributions there is good agreement between the data and simulated 
continuum events. In addition, by analyzing data taken on the T(4S) resonance, we observe 
the closely related semileptonic B decay with yields that agree well with predictions from 
the detector simulation. 

We determine the efficiency for detecting a supersymmetric B meson using Monte Carlo 
simulation. A range of values for the B mass, M(B), from 3.5 to 4.5 GeV/c 2 , has been 
explored. Also, the sneutrino mass, M(z>), was varied within the kinematically allowed range. 
In the process e + e~ — > bb, followed by hadronization, the supersymmetric B mesons have 
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energies E(B) distributed from M(B) to -Ebeam according to some unknown fragmentation 
function, which we approximate by a delta function 5(E(B) — Eq), with Eq given a value 
between M(B) and -Ebeam- We simulate both B mesons and their daughters, as well as 
the additional hadrons that result from fragmentation. We determine the dependence of 
efficiency on M(B), M(u), and Eq. To simulate the decay B — > D£~X, we have taken X to 
be a single pion and used three-body phase space for the decay. For b — > c£~v, we have used 
D£~u three-body phase space. 

These variations affect the B meson detection efficiency primarily through the lepton 
momentum spectrum. Table [I] lists a representative set of detection efficiencies and the 
resultant 95% confidence level upper limits on the cross section, assuming M(B) = 4.0 
GeV/c 2 , M{y) = GeV/c 2 , and E = E heam . The efficiency exhibits weak dependence on 
M(B), but lower values of Eq tend to soften the lepton momentum spectrum, resulting in 
reduced efficiency. For the smallest value of Eq we considered (3.7 GeV), the efficiencies are 
one-fourth of those with Eq = -E"beam- We give the upper limit on the product of the BB 
production cross section and the B — ► D^£{n, z>} branching fraction as a function of M(B) 
and Eq, assuming M(y) = GeV/c 2 . For all modes, the variation of the efficiency with 
M(p) is roughly linear between M(u) = GeV/c 2 and M(u) max = 0.6M (B) - 0.8 GeV/c 2 , 
where the efficiency vanishes. Hence, the upper limits for a finite M{y) would be those for 

M{v) = GeV/c 2 , scaled by a factor (l - -^g-)" 1 . 

The detection efficiency is also affected by the additional particles produced in fragmen- 
tation. From Monte Carlo studies, we have determined that for each fragmentation particle, 
the fractional decrease in efficiency is approximately 3% and depends on the momentum and 
species of the fragmentation particles. By considering Monte Carlo simulation of e + e~ — > cc, 
we estimate roughly two fragmentation particles per GeV of fragmentation energy. The 
detection efficiency has been corrected accordingly. A systematic error per fragmentation 
particle of ±2.5% for the D modes and ±5% for the D* modes has been included in the 
upper limits to allow for the uncertainty in this correction. 

In addition, the upper limits have been inflated to account for uncertainties in the detec- 
tor simulation and in the D and D* branching fractions, which amount to systematic errors 
of 20%. We also assign an error of 10% to the expected Standard Model yields due to under- 
tainties in the simulation of both fake and real leptons. The total systematic uncertainties 
are small compared to the statistical errors on the background-subtracted yields. 

Figure |3] shows the 95% confidence level upper limit on the product of a(e + e~ — ► BB) 
and B(B -> D^k{ir,v}) function of M(B) and Eq, assuming M(u) = GeV/c 2 . 

The strong rise of the upper limits with decreasing Eq/E^^ reflects the softening of the 
lepton spectrum. Because of the /3 3 dependence of the cross section near threshold, the cross 
section at \/s = 10.52 GeV is heavily dependent on the assumed M(B). These predicted 
cross sections are shown as the thick curves in Figure |3|. For both the Z)-dilepton and D*- 
dilepton signatures, the upper limits in nearly all of the kinematically allowed M(B)-E Q 
parameter space with M(P) = GeV/c 2 fall below the predicted bb production cross section. 
Below M(B) = 3.9 GeV/c 2 , the _D*-dilepton efficiency vanishes, so we exclude a B that 
decays B — > D*£{tt, v} only in the mass range 3.9-4.5 GeV/c 2 . No portion of the parameter 
space is excluded if M(u) is greater than 1.2 GeV/c 2 and 1.3 GeV/c 2 for the D-dilepton and 
.D*-dilepton signatures, respectively. 
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If the i?-parity-violating decay b — > c£~ were to occur, then some fraction of the time it 
would appear as the dressed decay B — > D£~ or B — > D*£~. We have examined the _D£~ and 
invariant mass distributions in D^*H + £~ events for evidence of a peak indicating such 
two-body decays. We find no evidence of a peak. Fitting the distributions to a polynomial 
background plus a Gaussian with width given by our experimental resolution (10 MeV/c 2 ), 
and stepping the Gaussian mean over the mass range, we obtain upper limits on the product 
of the bb production cross section and the B — > D^*'£~ decay branching fraction, as a function 
of M(B). These upper limits, which include 15% systematic errors on the yields, are shown 
in Figure [| For B — > D£~, the upper limit is less than 0.3 pb for all masses except near 
4.36 GeV/c 2 , where it is 1.0 pb. For B — > D*£~ the upper limit is weaker, typically 3 pb for 
masses around 4 GeV/c 2 , dropping to 0.3 pb for masses near 4.5 GeV/c 2 . 

In conclusion, we have searched for associated BB production in e + e~ collisions with 
center-of-mass energy below the BB production threshold. We assume a branching fraction 
of 100% for B — > D^£~n or B — > D^*H~u. Considering D-dilepton and D*-dilepton com- 
binations, we find no evidence of a light scalar bottom quark produced at y/s = 10.52 GeV. 
Upper limits on the BB production cross section depend on the assumed mass and energy 
of the B meson, as well as the mass of the sneutrino. For M{y) less than 0(1 GeV/c 2 ), 
the existence of a light scalar bottom quark with mass between 3.5 GeV/c 2 and 4.5 GeV/c 2 
has been excluded at the 95% confidence level. A light scalar quark decaying 100% of the 
time b — > c£i> and/or b — > c£ would have escaped our notice only if its decay matrix element 
results in a lepton spectrum much softer than three-body phase space. 
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TABLE I. Fitted D and D* yields for the .t)(*)-dilepton signature. The expected yields due to 
Standard Model processes are obtained from an analysis of simulated continuum events and include 
a systematic error of 10% due to uncertainties in the modeling of both fake and real leptons. The 
representative signal efficiencies given include the Z?M branching fractions and assume M(B) = 4.0 
GeV/c 2 , Eq = -Ebeam, and M{y) =0 GeV/c 2 . The corresponding 95% confidence level upper limits 
on the product of the e + e~ —* bb production cross section and the B —* D^£{ir,i)} branching 
fraction are calculated from the excess in the measured yield over the expected yield. These upper 
limits include systematic errors of 20% on the reconstruction efficiency. 



Mode 


Yield 


Expected Yield 


e S i 9 (%) <r(e + e" 


" -> bb) x B(B -» D^£{tt, v}) 


D°£+£- 


47.6 ± 20.0 


33.7 ± 14.7 ±3.4 


0.44 ± 0.02 




D ± £+£- 


37.4 ±25.1 


58.3 ±21.3 ±5.8 


0.86 ± 0.04 




D£ + £- 








< 2.7 pb at 95% C.L. 


D*°£ + £- 
D* ± (D°Tr ± )£+£~ 
D* ± iD ± 7r°)£+£- 

D*£+£~ 


4.9 ±3.4 
11.3 ±3.8 
0.2 ±3.6 


4.2 ±2.6 ±0.4 
4.0 ±2.0 ±0.4 
3.6 ±2.3 ±0.4 


0.08 ±0.01 
0.12 ±0.01 
0.09 ±0.01 


< 3.7 pb at 95% C.L. 
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FIG. 1. Invariant mass distributions for (a) D° — ► K~7r + and (b) D + — > A'~7r + 7r + candidates 
paired with a single lepton of either charge (open histogram) or with two oppositely charged leptons 
(shaded histogram) . The two -D-dilepton distributions have been scaled by a factor of 25 to facilitate 
comparison with the normalizing distributions. The fits of each distribution to a Gaussian and a 
linear background are shown in the dashed curves. 
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FIG. 2. Distributions of M(D*) - M(D) for (a) D*° -> D°ir°, (b) D*+ D°ir + , and (c) 
D* + — > D + tt° candidates paired with a single lepton of either charge (open histogram) or with 
two oppositely charged leptons (shaded histogram). The three Z)*-dilepton distributions have been 
scaled by a factor of 15 to facilitate comparison with the normalizing distributions. The fits of each 
distribution to a Gaussian and a quadratic background are shown in the dashed curves. 
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FIG. 3. Upper limits at the 95% confidence level on the product of the e + e~ — > BB pro- 
duction cross section and the B — > D^l~{ir, 9} branching fraction. The upper limits for the (a) 
-D-dilepton and (b) I?*-dilepton signatures are given for the full kinematic range of E(B) / -Ebeam 
as a function of the B mass, M(B), assuming M(y) = GeV/c 2 . The thick curves show the 
theoretically predicted cross sections, derived from the cross section for pointlike fermions and 

scaled by (3 s = (l - M{B) 2 /El ear ^j . For finite M(P), the upper limits are increased by a factor 



( 



M(v) 



The shaded histograms represent the largest upper limits, for each 



^ 0.6M(B)-0.8 GcV/c 2 

given M(B), on an i?-parity- violating B that decays B — > D^*>1 
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